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the  effects  of  inertial  forces  acting  on  the  CNS  through  the  freely  moving 
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SUMMARY  PAGE 
PROBLEM 


In  order  to  test  and  evaluate  Impact  protection  devices,  an 
impact-injury  model  for  restrained  humans  in  a  crash  environment  must 
be  developed.  Disruption  of  the  functioning  of  the  central  nervous 
system  (CNS)  is  an  important  consequence  of  impact  injury  involving 
the  head  and  neck,  and  is  an  important  consideration  in  the 
development  of  a  useful  impact-injury  model.  Ultimately, 
neurophysiological  criteria  for  functional  injury  of  the  CNS  are 
desired.  The  main  purpose  of  the  experiments  reported  here  is  to 
identify  some  of  the  measures  of  CNS  function  which  may  provide  the 
basis  for  establishing  such  criteria. 


FINDINGS 


The  experiments  reported  here  are  part  of  an  on-going  program  to 
test  the  neurophysiological  effects  of  indirect  or  inertial  head-neck 
impact  acceleration.  In  these  experiments,  eight  unanesthatized 
Rhesus  monkeys,  with  torsos  restrained  in  a  seated  position,  and  with 
head  and  neck  free  to  move,  were  subjected  to  peak  sled  accelerations 
in  the  -X  direction  ranging  from  42  m/s2  to  963  m/s2.  Recordings  of 
somatosensory  evoked  potentials  (SEP's)  were  made  using  recording 
electrodes  chronically  implanted  over  the  somatosensory  cortex,  and 
over  the  cervico-medullary  junction.  Electrical  pulse  stimuli  were 
delivered  at  a  rate  of  5  Hz  through  spinal  electrodes  located  at 
LI-L2.  SEP's  were  recorded  prior  to  impact,  through  the  Impact  event, 
and  subsequent  to  impact.  Qualitative  analyses  were  performed  on  both 
cervical  and  cortical  SEP's,  and  extensive  quantitative  analysis  was 
performed  on  cervical  SEP's.  The  results  of  these  analyses  indicate 
that  neurophysiological  indices  of  injury  may  include:  increases  in 
latencies  of  the  cervical  SEP  peaks  exceeding  2.5%;  large  changes  in 
the  amplitude  of  the  cervical  SEP;  changes  in  ripples  on  the  cortical 
primary  SEP;  and  substantial  and  persistent  changes  in  the 
surf ace-positive  cortical  primary  SEP.  In  particular,  analysis  of 
shifts  in  latency  of  the  cervical  SEP  suggests  the  possibility  of  an 
injury  threshold  in  the  vicinity  of  700  -  800  m/s2.  Smaller  shifts  in 
latency  occuring  near  600  m/s2  may  indicate  a  pre-injury  condition. 


RECOMMENDATIONS 


Quantitative  measures  of  neurophysiological  function  are  useful 
in  assessing  the  effects  of  inertial  forces  acting  on  the  CNS  through 
the  freely  moving  head  and  neck.  Continued  investigation  and 
extension  of  these  measurement  techniques  is  recommended. 
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I.  INTRODUCTION 


Impact  injury  involving  the  head  and  neck  disrupts  the  normal 
functioning  of  the  central  nervous  system  (CNS)  to  an  extent  dependent 
upon  the  severity  and  nature  of  the  trauma.  The  often  temporary  CNS 
dysfunction,  which  results  from  head  and  neck  trauma  (identified 
clinically  as  "concussion"),  is  of  special  interest  in  the  development 
of  a  useful  injury  model .  The  pioneering  work  of  Denny-Brown  and 
Russell  (5)  was  the  first  thorough  attempt  to  identify  the 
physiological  concomitants  of  direct  head  injury  and  led  to  numerous 
subsequent  Investigations  (e  ,  g . ,  7,  12,  13,  31).  These  studies 
shared  the  following  features:  (a)  the  Impact  blow  was  delivered 
directly  to  the  head  which  was  free  to  move,  (b)  anesthetized  animals 
(monkey,  dog,  cat)  were  subjects,  (c)  basic  vital  functions  and  EEG 
were  monitored.  Generally,  the  results  followed  a  typical  pattern, 
similar  to  the  one  described  by  Denny-Brown  and  Russell  (5)  and 
Williams  and  Denny-Brown  (33).  Severe,  but  non-fatal,  blows  resulted 
in  a  loss  of  corneal  reflex,  rise  in  blood  pressure,  fall  in  heart 
rate,  drop  in  EEG  amplitude  and  frequency,  sometimes  followed  by 
development  of  slow  waves.  These  effects  could  occur  in  the  absence 
of  any  apparent  brain  pathology  and  have  been  well  reviewed  (6,  21, 
30). 


In  the  decade  following  these  reports,  two  important  extensions 
to  these  basic  findings  were  made.  Foltz  and  Schmidt  (8),  using 
unanesthetized  monkeys,  stimulated  the  sciatic  nerve  and  recorded 
evoked  potentials  (EP's)  from  the  reticular  formation  (RF)  and  the 
medial  lemniscus.  In  six  out  of  eight  monkeys  receiving  severe  direct 
head  Impacts,  the  lemniscal  response  persisted  vdiile  the  RF  response 
was  abolished.  This  first  use  of  the  EP  in  a  head  injury  study 
indicated  that  the  non-specific  brainstem  gray  matter  was  functionally 
more  sensitive  to  impact  than  the  sensory  specific  ascending  pathway. 
Subsequently,  Friede  (9)  demonstrated  that  in  the  cat,  both  cervical 
stretch  and  a  blow  to  the  head  produced  the  same  loss  of  reflexes  as 
well  as  the  same  neuropathology  at  the  level  of  the  spinal  cord. 
This,  again,  suggested  that  the  lower  brainstem  might  be  the 
vulnerable  site  in  CNS  impact  dysfunction. 

More  recently,  Ommaya  and  his  co-workers  (10,  15,  18)  subjected 
both  the  unanesthetized  monkey  and  chimpanzee  to  non-impact  head 
acceleration  while  stimulating  the  median  nerve  and  recording  the 
somatosensory  EP's  at  the  cortex.  The  EP  amplitude  was  more  sensitive 
to  head  acceleration  than  the  EEG.  This  occured  in  animals  that  were 
rendered  unconscious  (loss  of  muscle  tone,  insensitivity  to  painful 
stimuli)  as  well  as  in  some  that  remained  conscious.  The  duration  and 
intensity  of  the  EP  effect  appeared  to  parallel  the  duration  of  the 
unconsciousness,  but  no  relationship  to  the  intensity  of  acceleration 
was  reported. 

These  earlier  results  provided  the  background  for  the  effort 
begun  in  late  1974  when  the  Naval  Biodynamics  Laboratory  (then  the 
Naval  Aerospace  Medical  Research  Laboratory,  Detachment  #1)  undertook 
the  first  of  many  experiments  designed  to  test  the  neurophysiological 
effects  of  indirect  or  inertial  head  acceleration.  In  these 
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experiments  the  restrained  torso  is  accelerated,  with  the  freely 
moving  head  and  neck  receiving  the  "indirect"  acceleration  through  the 
skeletal  and  soft  tissue  anatomy.  Results  from  the  early  experiments 
(2,  32)  indicated  that  the  cortical  somatosensory  EP  showed  a 
decrement  in  amplitude  and  an  increase  in  latency  following  non-lethal 
Impact.  These  changes  appeared  to  be  greater  with  increased  peak 
acceleration.  This  confirmed  the  utility  of  neurophysiological 
measures  in  assessing  the  effects  of  inertial  forces  on  the 
functioning  of  the  brain.  Ultimately,  neurophysiological  criteria  for 
functional  injury  to  the  CNS  are  desired.  The  main  purpose  of  the 
experiments  reported  here  is  to  identify  those  measures  of  CNS 
function  which  can  be  the  basis  for  establishing  such  criteria. 


II.  METHODS 


I I- A.  Experimental  Procedure 

The  adult  Rhesus  (Macacca  Mulatta,  ca.  10  Kg)  was  selected  as  the 
animal  model  and  the  somatosensory  system  as  the  pathway  for  testing 
CNS  function.  All  the  experiments  were  designed  to  use  unanesthetized 
animals,  restrained  in  a  sitting,  upright  position,  with  head  and  neck 
freely  moveable.  The  restraint  system  consisted  of  a  nylon  suit  which 
covered  the  entire  body  except  for  the  head  and  neck.  Straps  sewn  to 
the  suit  firmly  restrained  the  subject  to  a  fiberglass  chair  which  was 
molded  to  the  shape  of  the  subject's  back.  The  subject  was  seated  on 
a  410  kg  sled  \A\ich  was  accelerated  by  a  HXGE  system  with  a  one  meter 
stroke.  Peak  sled  accelerations  ranged  from  42  to  963  m/s2.  The 
subject  was  oriented  so  that  he  was  accelerated  in  the  -X  (24,  25) 
direction.  The  sled  was  decelerated  by  friction  (2  to  3  m/s2)  over  a 
distance  of  up  to  213  meters.  Precise  inertial  data  were  obtained 
from  an  array  of  transducers  rigidly  mounted  to  the  monkey's  skull. 
Physiological  data  included  EKG,  respiration,  and  cervical  and 
cortical  EP  activity;  the  exact  configuration  occasionally  varied. 
Average  heart  rate  was  computed  by  hand  from  EKG  strip  chart  records. 

The  eight  monkeys  used  as  subjects  for  these  22  experiments  were 
prepared  surgically  with  chronic  in-dwelling  electrodes  using 
procedures  previously  described  (28,  29).  Briefly,  strip  disc 
electrodes  were  implanted  epidurally  over  the  dorsal  spinal  cord  at 
-  L 2  and  over  the  cervico-medullary  (CV)  junction.  Bilaterally, 
subdural  electrodes  were  implanted  over  the  primary  somatosensory 
cortex  (SX) .  Three  months  were  allowed  for  recovery  from  implantation 
surgery . 

Constant  current  monophasic  pulses  at  a  nominal  rate  of  5  Hz  were 
delivered  to  one  pair  of  lumbar  electrodes.  The  stimulus  intensity 
ranged  from  0.3  to  3.0  ma  among  animals  and  was  selected  for  each 
animal  as  the  highest  intensity  consistent  with  the  apparent  comfort 
of  the  subject.  Large,  but  sub-maximal  EP's  were  recorded  from  the  CV 
and  SX  electrodes  through  a  telemetering  system  with  an  overall 
bandpass  of  30  to  1500  Hz  for  CV  and  10  to  1500  Hz  for  SX.  Table  1 
details  the  parameters  used  for  each  experiment. 


TABLE  1  SUMMARY  OF  EXPERIMENTAL  PARAMETERS 
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EP  data  aquisltlon  was  Initiated  30  to  45  minutes  prior  to 
Impact,  was  continued  through  Impact  delivery,  and  was  terminated  45 
to  90  minutes  after  Impact.  There  was  a  five  to  ten  minute  gap  In 
data  aquisltlon  ending  ten  minutes  prior  to  Impact.  The  data  were 
amplfied  on  the  sled,  telemetered  to  nearby  equipment  and  recorded  on. 
FM  tape.  The  data  were  digitized  off-line  on  a  hybrid  EAI  PACER  60(r^ 
computer  at  a  sampling  rate  of  at  least  20  kHz.  A  software-hardware 
design  was  used  vrtiich  synchronized  digitization  with  the  stimuli.  The 
digitized  data  were  then  processed  on  a  UNIVAC  1100^  series  computer. 
Some  of  these  data  were  also  analyzed  elsewhere  (19). 


II-B.  Data  Analysis 

Qualitative  analysis  was  performed  on  both  cervical  and  cortical 
EP's.  Starting  at  various  times  in  relation  to  impact,  average  evoked 
potentials  (AEP's)  with  N  ■  50  were  computed  from  sequences  of  EP's. 
These  AEP's  were  visually  examined  for  effects  related  to  impact. 

Quantitative  analysis  was  performed  on  cervical  EP's  only,  and 
was  used  to  determine  the  extent  to  iriiich  the  Impact  produced  shifts 
in  latency  of  various  peaks  of  the  cervical  AEP's.  Fig.  1  is  a  flow 
diagram  of  the  quantitative  data  analysis  procedure.  Two  types  of 
AEP's  were  computed:  "test"  and  "baseline".  Test  AEP's  were  each 
computed  with  N  -  10,  and  since  the  inter-stimulus  interval  was  0.2 
seconds,  each  spanned  an  interval  of  two  seconds.  For  each  run,  a 
sequence  of  test  AEP's  was  computed  starting  two  minutes  pre- impact, 
and  e  _ing  five  minutes  post- impact.  A  single  baseline  AEP  was 
computed  using  the  EP's  occuring  from  two  minutes  pre- impact  until 
impact.  N  was  approximately  580. 

Most  of  the  cervical  AEP  peaks  with  latencies  between  1.7  and  5.5 
ms  were  subject  to  latency  analysis.  Peaks  which  proved  too  unstable 
for  analysis  were  excluded.  In  order  to  determine  the  pre- impact 
latency  of  each  peak,  a  simple  peak  detection  algorithm  was  used.  For 
each  peak  of  the  pre-impact  test  AEP,  a  latency  interval  was  selected 
within  which  the  peak  appeared.  The  absolute  maximum  or  minimum  in 
this  interval  was  determined  and  used  as  the  peak  latency.  The  result 
was  discarded  for  that  test  AEP  if  the  absolute  extremum  was  at  the 
beginning  or  end  of  the  selected  interval.  For  each  peak,  the  median 
latency  for  all  pre-impact  AEP's  was  then  computed  and  used  as  a 
baseline  latency. 

To  determine  shifts  in  latency  produced  by  impact,  an  algorithm 
to  compute  normalized  cross-correlation  functions  (NCCF's)  was  used. 
Computation  of  an  NCCF  may  be  thought  of  as  a  search  for  a  selected 
portion  of  the  baseline  AEP  in  each  of  the  test  waves.  This  selected 
portion  of  the  baseline  AEP  is  called  the  "reference  wave"  and  was 
that  portion  of  the  baseline  AEP  between  two  specified  latency  limits 
on  either  side  of  the  peak  of  interest.  To  compute  the  NCCF,  the 
reference  wave  was  shifted  along  the  test  wave  am.  Pearson's 
correlation  coefficient  was  computed  between  the  reference  and  test 
wave  at  each  value  of  the  shift.  The  NCCF  is  the  correlation 
coefficient  as  a  function  of  the  amount  of  shift.  To  obtain  the  shift 
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FIGURE  1 

Flow  diagram  of  the  analysis  procedure  applied  to  cervical  EP  data 
from  each  experiment.  There  are  two  inputs:  1  -  one  large  block  of 
pre-impact  EP's;  2  -  numerous  3mall  blocks  of  EE's  from  both  pre-  and 
post- impact  periods.  There  are  four  outputs:  1  -  median  duration  of 
latency  shift;  2  -  median  latency  shift;  3  -  median  fractional  latency 
shift;  4  -  two-term  exponential  fits.  The  first  three  outputs  are 
single  values  plotted  in  fig.  4.  The  last  output  consists  of  graphs 
illustrated  in  fig.  5,  and  regression  coefficients  given  in  table  3. 
The  coefficients  for  one-term  exponential  regression  are  intermediate 
values  in  the  flow  diagram,  and  are  given  in  table  2. 
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in  latency,  each  NCCF  was  scanned  for  the  highest  value  of  the 
correlation  coefficient.  The  shift  in  latency  (of  reference  wave  with 
respect  to  test  wave)  at  this  maximum  was  taken  as  the  shift  in 
latency  for  the  test  AEP.  If  the  maximum  correlation  coefficient 
occurred  at  either  end  of  the  range  of  shift  in  latency  the 
corresponding  data  were  discarded.  The  NCCF  was  applied  to  the 
pre- impact  and  post- impact  test  AEP's. 

In  both  simple  peak  detection  and  NCCF  computation,  it  was 
necessary  to  select  a  latency  interval  over  which  the  AEP  was 
searched.  Selection  of  this  interval  was  based  on  both  physiological 
and  technical  considerations.  The  interval  was  made  long  enough  so 
that  shifts  in  latency  known  to  be  possible  from  extensive  visual 
examination  of  AEP  data  would  be  detected.  It  was  necessary,  however, 
to  keep  the  interval  short  enough  to  avoid  confusion  with  adjacent 
peaks.  Occasionally,  these  two  requirements  proved  incompatible,  in 
which  case  the  peak  was  not  used . 

Both  the  simple  peak  detection  and  the  NCCF  procedures  involved 
detection  of  an  extremum  over  a  sequence  of  points.  In  some  cases,  to 
increase  the  resolution  of  the  measurements,  quadratic  interpolation 
was  used.  The  best  parabola  (in  the  least-squares  sense)  was 
determined  for  an  odd  number  of  points  centered  at  the  previously 
determined  extremum.  The  extremum  of  this  parabola  was  then  taken  as 
the  final  result. 

Shifts  in  latency  for  each  peak  were  plotted  as  a  function  of 
time  over  the  experiment,  relative  to  impact  (figs.  2,  3,  5).  To 
compare  shifts  in  latency,  the  magnitude  of  the  effect  of  impact  on 
the  latency  of  each  peak  was  represented  by  two  variables.  The  first 
is  the  amplitude  of  the  effect,  i.e . ,  a  measure  of  the  deviation  of 
the  shift  in  latency  from  its  pre-impact  value  (not  to  be  confused 
with  the  amplitude  of  the  evoked  potential) .  The  second  is  the 
duration  of  the  effect. 

To  estimate  these  two  variables,  first,  values  of  shift  in 
latency  that  resulted  from  failure  of  the  NCCF  to  detect  the  selected 
peak  were  manually  removed.  Next,  amplitude  and  duration  were 
determined  by  regression  of  a  one  or  two  term  exponential  decay 
function  on  data  from  the  first  five  post-impact  minutes.  In  most 
instances,  the  initial  effect  of  Impact  was  an  increase  in  latency, 
and  in  about  85%  of  these  cases,  single  exponential  regression  was 
used.  However,  in  the  remaining  15%,  the  latency  shifted  to  values 
well  below  the  pre- impact  baseline  for  a  substantial  portion  of  the 
first  five  post-impact  minutes.  In  these  cases,  a  large  positive 
constant  was  added  to  the  post-impact  sequence  of  shifts  in  latency, 
and  two  term  exponential  regression  was  used.  This  generally  resulted 
in  one  term  with  a  time  coefficient  that  was  much  longer  than  the  five 
minute  analysis  interval,  and  could  thus  be  considered  essentially 
constant.  The  amplitude  of  the  effect  was  then  computed  as  the  sura  of 
the  two  exponential  amplitude  coefficients  and  minus  constant  added  to 
the  data.  This  final  amplitude  represented  the  initial  post- impact 
displacement  of  the  latency  from  the  pre-impact  baseline,  and  was,  in 
this  sense,  comparable  to  the  amplitudes  from  the  one  term  analyses. 
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The  shorter  of  the  two  time  coefficients  was  used  to  estimate 
duration,  although  it  was  not  always  comparable  with  those  from  the 
one  term  procedure. 

Aside  from  determination  of  the  amplitude  and  duration,  one  and 
two  term  exponential  regressions  were  also  used  to  determine  whether 
the  corresponding  decay  functions  wauld  provide  a  good  representation 
of  the  time  course  of  effects  of  Impact.  Only  data  from  high  level 
Impacts  were  used  for  this  purpose. 

For  exponential  regression  the  method  employed  was  a  modification 
of  the  modified  Gauss-Newton,  least-squares  method  recommended  by 
Metzler,  et  al .  (17,  pp.  3-9).  Initial  estimates  of  the  exponential 
coefficients  were  usually  obtained  with  a  modified  version  of  the 
program  CSTRIP  (20),  but  in  some  cases  initial  estimates  were  manually 
supplied. 

In  determining  the  least  squares  fit,  the  mean  square  deviation 
is  expressed  as  a  function  of  the  exponential  coefficients,  and  the 
absolute  minimum  of  the  mean  square  deviation  is  sought.  The  method 
used  to  locate  the  absolute  minimum  can  sometimes  be  "trapped"  in  a 
relative  minimum,  and  return  an  incorrect  result.  Generally,  this  was 
obvious  from  the  plots.  The  situation  can  often  be  corrected  by 
altering  the  starting  values,  or  dropping  a  few  data  points.  These 
methods  were  applied  where  appropriate.  In  other  cases,  the  correct 
least  squares  solution  did  not  involve  only  exponential  decay. 

Instead,  an  exponential  growth  term  with  a  short  time  coefficient 
appeared.  Such  results  were  excluded  from  consideration.  Detailed 
consideration  of  some  aspects  of  the  data  analysis  appears  in 
Appendix  B. 


III.  RESULTS 

Pre- impact  and  post- Impact  AEP's  from  all  experiments  analyzed 
appear  in  Appendix  A,  and  are  specified  by  experiment  number.  Each 
experiment  number  begins  with  the  characters  "LX".  Constant  latency 
cursors  are  provided  to  aid  in  the  observation  of  shifts  in  latency. 
Unless  otherwise  indicated,  the  results  reported  here  refer  to  data 
illustrated  in  Appendix  A. 

III-A.  Cortical  AEP's  —  Description 

The  cortical  recordings  were  bipolar.  However,  in  many  of  the 
subjects  the  recording  electrodes  were  selected  so  that  the  positive 
electrode  was  over  the  most  active  cortical  region  available,  while 
the  negative  electrode  was  relatively  inactive.  It  is  therefore 
reasonable  to  speak  of  surface  polarity  of  the  signal  as  if  the 
recordings  were  monopolar.  The  activity  observed  in  the  cortical 
recordings  may  be  divided  into  four  general  categories: 

1  -  Cortical  Slow  Waves.  A  positive  wave  began  at  about  5  to  12  ms, 
peaked  at  about  9  to  11  ms,  and  was  followed  by  a  negative  wave 
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that  peaked  at  about  15  to  24  ms.  The  positive  wave  was 
considered  to  be  the  first  cortical  response  to  the  afferent 
volley.  There  was  a  substantial  subject-to-subject  variation  in 
the  amplitude  of  the  cortical  slow  wave. 

2  -  Fast  Ripples.  These  appeared  on  the  rise  and  peak  of  the 

positive  slow  wave.  They  have  been  observed  in  numerous 
subjects  in  this  and  other  studies  after  stimulation  of  either 
the  dorsal  columns  or  the  median  nerve  of  Rhesus.  They  appear 
to  have  a  narrower  spatial  distribution  than  the  cortical  slow 
wave  and  under  some  conditions  the  cortical  slow  wave  can  change 
without  apparent  change  in  the  ripples.  These  ripples  may 
represent  afferent  activity  in  the  thalamo-cortical  radiations. 

3  -  Early  Fast  Activity.  This  appeared  at  practically  any  latency 

from  about  1  ms.  to  the  onset  of  the  cortical  slow  wave,  and  may 
have  been  due,  in  part,  to  far  field  activity.  In  some 
subjects,  an  especially  prominent  and  very  stable  diphasic  spike 
was  seen  with  about  3.4  ms  latency  after  spinal  stimulation.  It 
appears  that  this  fast  spike  is  localized  in  the  region  of  the 
primary  somatosensory  cortex.  It  may  be  an  antidromic  pyramidal 
tract  response. 

4  -  Early  Slow  Wave.  This  was  seen  in  a  few  of  the  cortical 

recordings  with  a  peak  latency  of  about  5  ms.  This  may  be  a  far 
field  response,  possibly  of  thalamic  origin. 


Ill— B.  Cervical  AEP's  —  Description 

The  cervical  recordings  were  bipolar,  and  no  significance  is 
assigned  to  the  polarity  of  the  waves.  The  earliest  peak  in  each 
subject  ranged  from  1.5  ms  to  2.9  ms.  Multiphaslc  fast  activity  of 
considerable  amplitude  typically  lasted  until  about  5  to  6  ms,  with 
many  low  amplitude  ripples  following.  Narrow  and  wide  "half-waves" 
can  be  described.  The  narrow  half-waves  were  roughly  .35  ms  wide. 

The  wide  half-waves  were  typically  two  to  three  times  as  wide  as  the 
narrow  half-waves,  and  often  had  either  a  double  peak,  or  an 
inflection  near  the  peak.  This  suggested  that  the  wide  half-waves 
consisted  of  two  nearly  overlapping  narrow  half-waves.  Where  a  double 
peak  was  observed,  high  impact  levels  sometimes  affected  the  two  peaks 
differently,  in  keeping  with  the  idea  that  two  "  components"  were 
present.  Each  "component"  may  represent  the  activity  of  fibers  with 
similar  conduction  velocities.  Variation  in  the  placement  of  the 
recording  electrodes  allows  each  "component"  to  appear  with  either 
polarity,  resulting  in  a  great  variety  of  wave  configurations. 


III-C.  Qualitative  Effects  of  Impact  —  Cervical  AEP's 

The  cervical  and  cortical  qualitative  results  were  derived  from 
visual  inspection  of  the  AEP's  illustrated  in  Appendix  A. 

In  five  experiments  on  three  subjects  (subject  AR0761, 


Somatosensory  Evoked  Potentials 


experiments  LX3009,  LX3010;  subject  AR8857,  experiment  1X3186;  subject 
AR0012*,  experiments  LX3693,  1X3694)  the  peak  sled  accelerations 
exceeded  700  m/s2  ("high  level  Impacts").  Cervical  EP's  were 
available  in  all  of  these  experiments.  In  two  subjects,  each 
receiving  two  high  level  impacts,  there  were  profound  changes  in  the 
cervical  AEP's.  In  one  subject  (subject  AR0761;  experiments  LX3009, 
LX3010)  the  cervical  AEP's  were  nearly  or  entirely  obliterated  by 
Impact.  At  796  m/s2  (1X3009)  the  AEP  waveshape  remained  distorted  for 
at  least  60  seconds.  At  963  m/s2  (LX3010)  the  pre- impact  waveshape 
was  certainly  not  restored  at  two  minutes,  and  probably  not  at  four 
minutes  post- Impact.  In  the  other  subject  (subject  AR0012, 
experiments  1X3693,  1X3694)  in  which  strong  effects  were  observed,  the 
early  portion  of  the  cervical  AEP  (2.11  ms  and  earlier)  was 
drastically  reduced  in  amplitude,  while  a  later  major  peak  (3.44  ms) 
did  not  change,  or  possibly  Increased  in  amplitude.  In  the  remaining 
experiment  (subject  AR8857,  experiment  1X3186,  810  m/s2)  there  is  no 
obvious  effect  in  the  cervical  AEP.  In  the  four  experiments  in  which 
strong  effects  were  observed,  there  were  obvious  increases  in  peak 
latencies  observable  as  the  AEP  components  became  recognizable  after 
Impact.  In  1X3186  visual  inspection  showed  no  increase  in  latency, 
although  quantitative  analysis  using  AEP's  with  an  N  of  10  (see  below) 
indicated  unequivocally  that  there  was  a  considerable  increase  in 
latency  of  short  duration. 

In  five  experiments  on  three  subjects  (subject  AR8824,  experiment 
LX3184;  subject  AR8802,  experiments  1X3705,  1X3706;  subject  ARNA02, 
experiments  LX3714,  LX3715)  the  peak  sled  acceleration  was 
approximately  600  m/s2  ("middle  level  impact").  Cervical  AEP's  were 
available  in  all  three  subjects,  although  in  one  subject  they  were 
marginal  (subject  AR8824,  experiment  1X3184).  There  were  no  profound 
effects  on  the  cervical  AEP  amplitude  although  some  small  amplitude 
effects  can  be  observed  (subject  AR8802,  experiments  LX3705,  1X3706). 
There  was  no  apparent  effect  on  waveshape  with  one  noteworthy 
exception.  This  was  the  early  portion  of  the  cervical  AEP  in 
experiment  LX3706,  and  will  be  discussed  below.  Some  Increases  in 
peak  latency  can  be  seen  in  experiments  1X3705  and  1X3706. 

The  peak  sled  accelerations  of  the  remaining  experiments  (see 
table  1)  were  below  440  m/s2  ("low  impact  levels").  Substantial 
changes  in  cervical  AEP's  did  not  occur  in  these.  There  are 
occasional  shifts  in  latency  of  the  cervical  AEP  peaks,  but  these  were 
small.  No  apparent  AEP  amplitude  effects  related  to  low  level  Impact 
were  noted  by  visual  examination  of  these  AEP's. 

I 


*  Histological  examination  revealed  old  traumatic  alterations  in  the 
cervico-medullary  region  of  subject  AR0012  (27).  The  results  from 
this  subject  should  therefore  be  considered  tentative. 
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CERVICAL  AEP  LATENCY  SHIFT 

RUN  LX 3009 >  796  M/S?  PEAK  3.45  MS.  N-50 


lit  put 


FIGURE  2 

Linear  and  three-term  exponential  regression  procedure  applied  to 
shifts  of  cervical  AEP  peaks  measured  by  normalized  cross-correlation 
for  all  of  experiment  LX3009.  AEP’s  were  computed  with  N  *  50,  and 
except  in  the  region  of  Impact,  only  every  sixteenth  EP  was  used.  For 
the  purpose  of  regression,  the  data  were  weighted  by  the  actual  time 
covered  by  the  AEP  (length  of  the  horizontal  line  representing  the 
data  point).  First,  linear  regression  was  applied  to  the  pre-impact 
data  resulting  in  a  slope  of  -35  ps/hr  and  a  Y  intercept  of  -12  ps. 
This  was  extrapolated  Into  the  post-impact  period  and  subtracted  from 
the  data.  Next,  the  constant  200  ps  was  added  to  the  post-impact 
data.  Finally,  three-term  exponential  regression  was  applied 
resulting  In  the  coefficients: 


AMPL  (ps): 

TIME  COEF  (min): 


551. 

-.356 


The  first  term  has  a  value  of  214  ps  at  Impact,  and  191  ps  at  40  min. 
post-impact.  Subtracting  the  200  ps  constant  originally  added  to  the 
data  to  the  first  term,  the  contribution  of  the  first  term  ranges  from 
+1 A  ps  at  Impact  to  -9  ps  at  40  min.  post-impact.  The  long  time 
coefficient  of  the  first  term  means  the  first  term  is  essentially 
constant  with  its  amplitude  range  being  within  the  range  of 
experimental  error.  This  term  accounts  for  the  drift  below  the  linear 
regression  line  at  the  end  of  the  experiment.  The  other  two  terms 
account  for  most  of  the  curvature.  The  final  weighted,  post- Impact 
RMS  deviation  was  10.9  ps. 


PEAK  3.17 


IMPACT 


0.20  MSEC 


20.  MIN 


PEAK  4.58 


FIGURE  3 


Plots  of  the  shifts  in  latency  at  two  peaks  measured  with  NCCF's  for 
all  of  experiment  LX3706.  AEP's  were  computed  with  N  -  100.  The 
increase  in  latency  of  the  peak  at  4.58  ms.  Induced  by  impact 
Illustrates  the  most  common  finding.  The  substantial  decrease  In 
latency  of  the  peak  at  3.17  ms.  is  unusual  and  was  caused  by  the 
Impact- Induced  appearance  of  a  new  component  that  overwhelmed  the 
original  peak. 
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III-D.  Qualitative  Effects  of  Impact  —  Cortical  AEP's 

Cortical  AEP's  were  available  in  two  of  the  three  subjects 
receiving  high  level  impacts  (subject  AR0761,  experiments  LX3009, 
LX3010;  subject  AR8857,  experiment  LX3186).  The  most  consistent 
finding  in  these  three  experiments  was  that  the  ripples  on  the  rise 
and  peak  of  the  positive  cortical  primary  response  were  reduced  by 
Impact.  It  appeared  that  the  amplitude  of  the  ripples  covaried  with 
the  amplitude  of  the  cervical  AEP.  The  cortical  primary  itself  was 
affected  in  varying  degrees  by  the  Impact,  corresponding  to  the 
strength  of  the  cervical  effect. 

Cortical  AEP's  were  obtained  from  two  of  the  subjects  receiving 
middle  level  impacts  (subject  AR8802,  experiments  LX3705,  LX3706; 
subject  ARNA02,  experiments  LX3714,  LX3715).  Distinct  ripples  were 
not  seen,  and  in  one  subject  (AR8802)  the  cortical  AEP  amplitude  was 
low.  A  substantial  decrease  in  the  amplitude  of  the  cortical  AEP 
occured  in  the  same  subject  in  which  cervical  Increases  in  latency 
were  seen  (subject  AR8802,  experiments  LX3705,  LX3706) . 

Substantial  changes  in  cortical  AEP's  were  not  seen  in  the  low 
level  experiments.  Cortical  ripples,  where  present,  were  never 
altered.  There  were  no  apparent  cortical  AEP  amplitude  effects 
related  to  low  level  Impact. 


III-E.  Quantitative  Analysis  of  Cervical  Shifts  in  Latency 

The  shifts  in  latency  of  several  cervical  EP  peaks  were  analyzed 
using  the  NCCF  and  exponential  regression  procedures  (see  Appendix  B 
for  details).  This  analysis  was  done  using  AEP's  with  10  EP's  each, 
and  spanning  a  seven  minute  time  Interval  starting  at  two  minutes 
pre-impact.  The  reason  for  this  restricted  analysis  Interval  is 
illustrated  in  figs.  2  &  3.  Fig.  2  illustrates  a  nearly  linear 
decrease  in  latency  over  the  entire  course  of  the  experiment, 
including  the  pre-impact  period.  Fig.  3,  on  the  other  hand, 
illustrates  shifts  in  latency  which  exhibited  an  eratic,  slow  drift 
during  the  experiment.  Such  slow  drifts  were  often  seen,  and  made  it 
possible  to  follow  reliably  the  impact- induced  shifts  in  latency  for 
no  more  than  about  ten  minutes  after  Impact. 


III-F.  Amplitude  and  Duration  of  Shift  in  Latency 

Table  2  summarizes  the  results  of  estimation  of  the  amplitude  and 
duration  of  shifts  in  latency  produced  by  Impact  for  all  available 
cervical  AEP  peaks  in  the  range  1.7  ms  to  5.5  ms.  The  shifts  were 
measured  using  the  NCCF.  Each  peak  is  identified  by  its  pre-impact 
median  latency  (LAT)  in  milliseconds  obtained  from  AEP's  computed 
using  data  collected  during  the  two  minutes  immediately  prior  to 
impact.  Alignment  of  comparable  peaks  into  columns  is  only 
approximate  across  subjects.  However,  within  subjects,  the  alignments 
are  correct,  and  pre-impact  latency  differences  within  subjects 
represent  actual  shifts  in  latency  between  experiments. 
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For  each  exponential  regression,  the  exponential  amplitude  (AMP) 
is  given  in  ps  of  increase  in  latency,  and  the  time  coefficient  (TIME) 
is  given  in  minutes.  A  negative  time  coefficient  indicates 
exponential  decay. 

For  plotted  data  which  had  no  visually  apparent  exponential 
patterns,  there  were  two  common  outcomes  of  the  exponential 
regression.  First,  a  long  time  coefficient  resulted  and  the 
exponential  amplitude  represented  an  estimate  of  the  mean  of  the  data. 
Second,  a  short  time  coefficient  resulted  and  the  exponential 
amplitude  represented  latency  shifts  in  the  early  post- impact  part  of 
the  data.  Both  types  of  results  appear  in  the  table.  In  each  case, 
the  time  coefficient  is  meaningless.  Such  results  were  Identified  by 
visual  examination  of  the  graphs,  and  the  corresponding  time 
coefficients  appear  in  parentheses.  Where  the  exponential  algorithm 
returned  no  useful  result,  and  visual  inspection  of  the  graphs 
revealed  no  apparent  effect  due  to  Impact,  zeros  were  filled  in. 

Two  points  should  be  emphasized  in  regard  to  this  table.  First, 
the  initial  few  post- impact  points  were  generally  discarded. 
Additionally,  in  high  level  Impact  experiments,  for  the  first  10  to  60 
seconds  post-impact,  the  shifts  in  latency  were  very  variable.  This 
resulted  in  a  substantial  number  of  shifts  in  latency  in  this  Interval 
being  considered  outliers,*  and  being  discarded.  The  exponential 
amplitudes  represent  extrapolation  of  the  shifts  in  latency  back  to 
impact  from  the  post- Impact  period,  and  do  not  represent  actual 
observations  of  the  shifts  in  latency.  The  actual  values  of  the 
maximum  shifts  in  latency  that  occured  as  a  result  of  Impact  are  not 
known  due  to  the  severe  distortion  of  the  evoked  potentials  as  well  as 
the  artifact  produced  by  Impact.  The  second  point  is  that,  the 
one-term  exponential  regression  is  used  as  a  good  means  of  providing 
objective  comparison  among  experiments,  and  not  as  a  model  for  the 
actual  underlying  process  (see  Appendix  B) . 

The  results  in  table  2  are  summer ized  in  fig.  4.  For  each 
experiment,  a  single  median  value  computed  from  all  of  the  peaks  for 
that  experiment  is  shown  as  a  point  in  each  of  three  plots.  The 
computations  for  each  plot  were  as  follows: 

4A  -  The  median  of  the  negative  time  coefficients  that  are 

meaningful  (not  in  parentheses)  and  are  not  associated  with 
negative  exponential  amplitudes. 

4B  -  The  median  of  non-negative  exponential  amplitudes. 


*  In  general,  an  "outlier"  is  defined  here  as  a  data  value  that 
deviates  enough  from  other  data  values  so  that  it  is  likely  to  contain 
erroneous  information.  Operationally,  the  definition  of  "outlier"  is 
embodied  in  the  procedure  that  is  used  to  exclude  it.  Such  procedures 
will  always  Imply  assumptions  about  the  nature  of  the  data,  (see 
Appendix  B) 


Somatosensory  Evoked  Potentials 


o  S 


sc  a  o 

<  h  * 

3  O  M  £ 

o  <  o  ft 

“  *  <  * 

a  »■  j  oo 


§  8 

*  *  5  - 

S  ff  £ 
a  3  3 


g  *  -6 
3  1 


c  a  > 

o  «i 

■H  £  H 

4->  O 

o  to  ® 

C  41  ^ 

9  o 

4*  C 

O  44 

w  a 

*j  w 
«  c  o 

«  *  8 
►.  o. 
o  *o 

b  «  « 

w  r-t  «-> 

«j  oo  a 

«  B  « 

H  >H  CO 

co  a) 

4-1  o 

0(0  0. 

0) 

to  s-.  o 

*•>  .e 

<44  4-1 

>rl  T>  C 
X  V  V 
■ 

B  ®  O 
O  TO  ® 

to  a 
cox 
o  a  o) 
■h  v 

TOO® 

w  x 

V  Q)  4J 

o  o 

60  TO  O 
®  O 

O  44  >44 

a!: 

o  o  .-c 
C  ®  X 

®  5.  TO 

C  X  4-> 

o  ® 

cl  c 

K  X  *4 

V  o 

«  TO 

E 

2  8  8. 

I  o 

«  <44  rl 
O  TO  *c3 

44 

grl  O 

2  ° 
O  TO  44 
4-4  0) 

o  ® 

T3  9 

V  ®  >-4 

h 

V  B 

TO  •  TO 
6  <H 
CD  O  TO 
TO  "H  TO 
9  0  0 
i-4  TO 
TOO® 
>  ®  X 

1-4  *4 
«  ® 

®  a  » 
o  u 
X  TO 


B 

44  TO  «  B 

o  4j  XI  o 

TO  44  Tt 

B  5  TO 
O  0  X  J  " 

■H  J3  ®  TO 
4-4  ®  O  4-4 

TO  X  M  OOr-t 
O  *4  TO  «  9 
TO  «  O  TO 
TO  I  ft  « 

o  o 

®  CJ  X  TO 
X  O  TO  TO 
*4  TO  B  TO 
•  TO  *.-<  TO 

op  y  1-4  x 

G  44  o  4-> 

t4  tI  O  » 

4-1  X  <44  O  X 

B  «  *4 

®  ►*’9  rj 

to  u  b  a 
®  o  a  • 

OB®  r* 

O.  ®  44  po  « 
V  4J  TO  *“«* 

O  TO  i-4  B  0 

f~l  M 

44  C  O 

AH'-*  O 

®  TO  ‘to 

i4  -rt  44  O 
u  44  44  TO  B 
■t4  "O  -H  -r4  TO 
44  6  £  TO  C 

44  "H  CO  U  44 

0)  0 

O  B  r-l  O 


m 

Q  9  Si 
ill 


O  1-4  O 
TO  S. 
TO 

*4  J£  X 

O  TO  o 
-4  ®  TO 
04  CL  Cd 


O  B  r4  O 
O  TO  TO  *9  TO 

t4  -r4  TO  44 

®  TO  44  44  TO 
0  TO  -H  44  « 

■3  0  a  o  o 

44  4  H  Ml 

«  a 

®  X  ®  TO 

>  44  £.  «  a 

•H  44  jC  o 
44  |  44  4-1 

8)  CO  o'  00  TO 

®  -h  a  o 
B  *9  ® 

■  TO  4  H 

a  »  >  a  « 
TO  0)  TO  a  o 

TO  O  TO  8  O 

-a  ®  o  « 

«  X  >s  O 

0  44  x>  TO3 
B  ®  V 

TO  TO  "9  O  r-4 

X  O  TO  O  « 

44  TO  B  04 
CL  t-C  TO  ft 

I  TO  X  TO 

<TO  X  J!  ?• 
O  TO 

TO  TOO 

•  CO  ►,  ex  o 

TOO®  <44 

44  x  a  « 

<44  44  TO  X  TO 

14  44  44  44 

JB  60  TO  B 

TO  0  r-5  O  *4 

■H  O  Q 
►.•9  flo  8. 
0  9  -rC 

B  i-t  ►»  TO 

V  O  44  O  TO 
44  X  <44  C  Q 

TO  «  iH  TO  .B 

<-4  £44 

6-.  TO  TO 
TOO  r-C  0 

►  B  _  44 

4-1  TO  TO  B 

4  44  B  TO  >1 

8,3 

TO  44  TO  o 

fi  B  O  0 

•rt  TO  *9 

£2  O  44  ® 

4  4  <4  U  4 

XTOsa. 

»  CO  <3  -3  TO 

J0  44  | 

TO  TO  4-)  TO  TO 

8.5  5  8.  S 


O  B  S 
TO 

a  o  o 

0CO 

•3  <r  m 


Somatosensory  Evoked  Potentials 


4C  -  The  median  of  non-negative  exponential  amplitudes  each 
divided  by  the  corresponding  pre- impact  median  latency. 

Fig.  4A  shows  that  the  median  time  coefficients  were  very  variable 
even  at  low  impact  levels,  and  there  was  little  if  any  relationship 
between  the  median  time  coefficient  and  Impact  level.  Figs.  4B  &  4C, 
on  the  other  hand,  show  that  there  was  a  strong  relationship  between 
the  exponential  amplitudes  (initial  shift  in  latency)  and  Impact 
level.  At  Impact  levels  below  about  500  m/s2,  there  were  small 
positive  median  latency  shifts  (mean  -  40  ps)  and  small  positive 
median  fractional  latency  shifts  (mean  -  9.3  ps/ms).  At  about  600 
m/s2  there  was  a  tendency  for  shifts  in  latency  following  Impact  to 
exceed  those  at  lower  levels.  Above  700  m/s2,  impact  produced 
increases  in  latency  that  exceeded  those  at  lower  levels.  There  are 
not  enough  data  to  determine  the  form  of  the  functional  relationship 
between  Impact  and  the  initial  shift  in  latency.  As  a  first 
approximation  to  the  curve  at  high  impact  levels,  linear  regression 
was  applied  to  the  results  in  figs.  4B  and  4C,  using  only  those  values 
with  impact  levels  above  500  m/s2.  The  results  were: 

figure  slope  impact  intercept 

4B  6.8  p8/(m/82)  490  m/s2 

4C  2.3  (ps/ms)/(m/s2)  500  m/s2 


II1-G.  Two-Term  Exponential  Fits. 

Only  in  the  high  impact  (727  m/s2  and  above)  experiments  were  the 
data  generally  noise  free  enough  to  allow  meaningful  comparisons 
between  one-term  and  two-term  regression.  The  results  are  presented 
in  table  3.  For  each  peak,  two  pairs  (exponential  amplitude  and  time 
coefficient)  of  results  are  shown  representing  each  of  the  two 
exponential  terms.  Also  shown  are  the  percent  decreases  in  mean 
square  deviation  obtained  with  the  two-term  regression  compared  with 
the  one-term  regression.  Substantial  Improvement  is  often  obtained. 
This  is  illustrated  in  fig.  5,  where  the  decrease  in  mean  square 
deviation  represents  a  clear  improvement  in  the  fit.  As  in  the  case 
of  the  single  exponential  results,  the  exponential  amplitudes 
represent  extrapolations. 


I1I-H.  Heart  Rate  Effects 

For  each  experiment  in  table  2,  the  average  heart  rate  was 
compared  for  two  20  second  intervals:  one  immediately  before,  and  one 
Immediately  after  Impact.  The  difference  between  these  two  rates 
("relative  bradycardia",  23)  is  plotted  in  fig.  6  as  a  function  of 
peak  sled  acceleration.  A  linear  regression  using  all  of  the  data 
points  is  also  plotted,  and  has  a  slope  of  -0.14  (beats/min)/(m/s2) . 
The  impact  Intercept  (zero  change  in  heart  rate)  is  113  m/s2.  The 
duration  of  the  measured  bradycardia  did  not  exceed  one  minute  except 
in  experiment  LX3184,  where  a  slight  relative  bradycardia  persisted 
for  approximately  2.5  min.  post- Impact. 
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Somatosensory  Evoked  Potentials 


CERVICAL  AEP  LATENCY  SHIFT 


POLY- EXPONENTIAL  REGRESSIONS 

(ON  132  OF  145  POST-IMPACT  POINTS) 


RUN  LX3009,  796  M/82,  PEAK  6.46  MS.  N-10 


nqjAh"^J>vr|L|f 


IMPACT 


260  U8KC 
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CERVICAL  AEP  LATENCV  SHIFT 


POLY-EXPONENTIAL  REGRESSIONS 

(ON  124  OP  151  POST-IMPACT  POINTS) 


RUN  LX3694, 


767  M/S2.  PEAK  1.77  M8.  N-10 
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,Af  va — r— j — in> 

FIGURE  5 

A  and  B  are  each  plots  of  the  shift  in  latency  for  one  peak  (indicated 
by  the  peak  latency  over  the  plot)  in  high  Impact  experiments.  Each 
is  from  a  different  subject.  AEP's  were  computed  with  N-10,  and  the 
interval  shown  ranges  from  2  min.  pre- impact  to  5  min.  post-impact. 
These  are  the  parameters  that  were  used  for  the  quantitative  analysis 
in  this  report.  The  shifts  for  each  peak  are  illustrated  twice, 
showing  one  and  two  term  exponental  regression  curves.  In  each  case 
the  two  term  fit  is  substantially  better  than  the  one  term  fit. 
Pre-impact  linear  regression  (Appendix  B)  was  not  used  in  these  fits. 
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POST-IMPACT  CHANGE  IN  HEART  RATE 


PEAK  SLED  ACCELERATION 
I  M/S  2I 


FIGURE  6 

Change  in  heart  rate  plotted  as  a  function  of  peak  sled  acceleration. 
Heart  rate  was  computed  from  strip  chart  EKG  records  for  the  20 
seconds  immediately  preceeding  Impact  and  for  the  20  seconds 
immediately  following  impact.  The  change  in  heart  rate  is  the 
post- Impact  rate  subtracted  from  the  pre-impact  rate.  The  average 
pre-impact  heart  rate  across  experiments  was  216  beats/min  with  a 
standard  deviation  of  27  beats/min.  A  linear  regression  is  also 
plotted  and  has  a  slope  of  -0.14  (beats/min) /(m/s2)  and  an  Impact 
intercept  of  113  m/s2. 


An  unusual  change  occurea  in  experiment  LX3706  (624  m/s2).  A 
substantial,  long-lasting  decrease  in  latency  was  seen.  The 
significance  of  this  result  is  seen  in  Appendix  A,  experiment  LX3706. 
At  the  time  of  impact,  an  early  EP  component  which  previously  had  been 
present  with  low  amplitudes,  if  at  all,  appeared  with  high  amplitude. 
This  additional  component  substantially  distorted  the  peak  at  3.17 
milliseconds,  resulting  in  the  apparent  decrease  in  the  latency  of  the 
peak.  The  impact  apparently  lowered  the  stimulus  threshhold  for  this 
component.  Fig.  3  (upper  trace)  presents  the  time  course  of  the 
recovery  of  this  effect  over  the  entire  experiment.  The  effect  begins 
to  decay  at  a  rate  consistent  with  the  -3.2  minute  time  coefficient 
given  in  table  2,  but  at  about  five  minutes  post-impact,  the  trend 
reversed  and  the  new  early  component  again  gained  strength. 
Subsequently  the  shift  in  latency  followed  the  variation  in  the  peak 
at  4.58  milliseconds  which  was  not  comparably  distorted.  In  the 
absence  of  such  peak  distortion,  substantial  decreases  in  latency  as  a 
consequence  of  impact  have  not  been  observed. 


IV.  DISCUSSION 

IV-A.  Amplitude  and  Duration  of  Shift  in  Latenc; 

The  one-term  exponential  regression  was  used  to  provide  objective 
measures  of  the  amplitude  and  duration  of  the  shifts  in  latency  of 
various  peaks  of  the  cervical  EP's  produced  by  impact.  The  data  show 
that  there  is  a  positive  relationship  between  the  strength  of  impact 
and  the  increase  in  latency  of  the  cervical  AEP.  The  data  are 
insufficient  to  determine  the  mathematical  form  of  the  relationship 
between  impact  strength  and  increase  in  latency.  Linear  regression 
was  used  as  a  first  approximation  to  the  high-impact  portion  of  the 
curve  (figs.  4B  and  4C) .  The  linear  regression  shows  that  there  was 
no  change  in  the  shift  in  latency  until  500  m/s2.  This  may  be 
contrasted  with  the  results  of  Weiss  and  Berger  (32)  who  found  that  a 
complex  measure  of  cortical  EP's  changed  at  all  levels  of  Impact. 
Further  analysis  of  the  cortical  EP  data  from  these  experiments  should 
provide  additional  information  in  this  regard. 

In  most  of  the  high  level  impact  experiments  (727  m/s2  and 
higher),  the  cervical  AEP  was  either  substantially  modified  or 
virtually  obliterated  (Appendix  A:  subject  AR0761,  experiments  LX3009, 
LX3010;  subject  AR00L2,  experiments  LX3693,  UC3694;  but  not  subject 
AR8857,  experiment  LX3186).  Under  such  conditions,  the  NCCF  algorithm 
operates  marginally,  and  the  variability  of  the  measurements  of  shifts 
in  latency  is  high.  For  the  purpose  of  exponential  regression,  many 
of  these  shifts  in  latency  were  considered  outliers,  and  were 
eliminated.  This  leaves  open  two  questions  of  considerable  interest. 
First,  what  is  the  greatest  latency  shift  that  can  actually  be 
produced  by  impact?  Second,  does  impact  Increase  the  variability  of 
the  latency?  Extensive  visual  examination  of  the  AEP's  showed  a  low 
signal-to-noise  ratio  in  much  of  the  AEP  data  immediately  following 
high  level  impact.  Therefore,  given  the  small  number  of  relevant 


Somatosensory  Evoked  Potentials 


experiments,  these  issues  cannot  be  resolved  at  present,  but  may  be 
after  additional  experimentation. 


IV-B.  Two -Term  Exponential  Fits 

The  double  exponential  model  provides  an  excellent  fit  to  the 
shifts  in  latency  for  the  first  five  post-impact  minutes  in  high  level 
Impact  experiments.  However,  numerous  other  models  with  four  or  fewer 
free  coefficients  are  likely  to  fit  these  results  at  least  as  well. 
This,  as  well  as  the  fact  that  the  double  exponential  results  rely  on 
only  four  experiments  on  two  subjects,  leaves  a  double  exponential 
hypothesis  tentative. 

This  hypothesis  is  that  high  level  impact  produces  two  different 
increases  in  cervical  EP  latency,  each  decaying  at  its  own  rate.  The 
two  hypothetical  effects  may  be  differently  related  to  impact.  It 
should  be  noted  that  despite  the  possible  inaccuracy  of  the  longer 
time  coefficients,  the  use  of  the  two-term  regression  results  in  an 
increase  in  accuracy  of  the  shorter  time  coefficient.  This  is  because 
of  the  better  fit  to  the  early  post-impact  part  of  the  curve  in  the 
two  term  fits.  The  shorter  time  coefficient  ranges  from  -14  to  -66 
seconds,  and  the  longer  time  coefficient  from  -102  seconds  to  very 
large  negative  values.  Assuming  the  double  exponential  hypothesis  is 
correct,  time  coefficient  magnitudes  up  to  about  ten  minutes  are 
reliable. 

Exponential  amplitudes  for  changes  in  latency  as  high  as  860  ps 
were  obtained  with  two-term  regression  in  the  high  level  impact 
experiments  (table  3).  However,  the  greatest  shift  in  latency  found 
by  direct  visual  examination  of  cervical  AEP's  in  these  experiments 
was  from  300  to  700  ps.  One  possible  explanation  for  this  discrepancy 
is  that  immediately  after  Impact  the  EP  disappeared,  and  after  some 
delay  reappeared  with  an  increased  latency  which  then  began  to  decay. 
Introduction  of  such  a  delay  into  the  regression  procedure  would 
reduce  the  exponential  amplitudes. 


IV-C.  Heart  Rate  Effects 


The  bradycardia  observed  after  impact  was  also  reported  by  Taylor 
and  Rhein  (23)  in  human  subjects  experiencing  indirect  head-neck 
acceleration.  At  a  peak  acceleration  of  15g  they  found  very  slight 
decrements  in  heart  rate  after  -X  acceleration,  and  substantial 
decrements  after  +X  acceleration.  The  results  of  our  experiments 
using  -X  acceleration  on  Rhesus  (fig.  6)  suggest  an  Inverse 
relationship  between  heart  rate  and  peak  sled  acceleration.  There  is 
a  lack  of  any  apparent  threshold,  in  contrast  to  the  cervical  AEP 
latency  results  (fig  4,  BSC).  The  extent  to  which  these 
cardiovascular  effects  are  related  to  the  measures  of  CNS  function 
remains  to  be  elucidated. 
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IV-D .  Possible  Temperature  Effects 

Stockard,  et  al .  (22)  report  the  relationship  between  the 
relative  latencies  of  peaks  of  the  auditory  far  field  potentials  and 
esophageal  temperature  in  human  patients  undergoing  anesthesia  and 
induced  hypothermia.  Using  their  data  from  37.1°C  and  34.5°C,  we 
computed  a  temperature  dependence  of  -38  (ps/ms)/°C  ((microseconds  of 
shift  per  millisecond  peak  latency)  per  degree  centigrade).  Marshall 
and  Donchin  (16)  report  the  relationship  between  circadian  variation 
of  oral  temperature  and  auditory  far  field  potentials  in  normal  human 
subjects.  From  their  figure  1,  and  tabular  data  on  "Latency  to  peaks" 
we  computed  a  temperature  dependence  of  -30±7  (ps/ms)/“C.  If  the 
cervical  EP's  were  to  exhibit  comparable  temperature  dependencies,  the 
slow  drifts  observed  (figs.  2  &  3)  could  be  due,  at  least  in  part,  to 
changes  in  core  temperature  during  the  experiment.  Further  work  is 
required  to  clarify  this  matter.  The  influence  of  core  temperature 
may  place  an  upper  limit  on  the  time  interval  over  which  the  latency 
data  can  reliably  be  used  to  determine  effects  of  impact.  However, 
measurement  of  the  sudden  effects  of  impact*  will  not  be  contaminated 
over  a  five  to  ten  minute  interval,  since  the  core  temperature  cannot 
change  rapidly  enough. 


V.  CONCLUSIONS 


Decreases  in  the  amplitude  of  cortical  EP's  have  been  observed  as 
a  result  of  direct  impact  (10,  15,  18)  and  indirect  impact  (3,  32). 
However,  there  is  evidence  that  large  amplitude  changes  resulting  from 
cognitive  factors  may  occur  in  primary  EP's  recorded  from  both  the 
cortex  and  the  lower  brainstem  (e.g. ,  l,  14).  While  shifts  in  EP 
latency  under  such  conditions  are  less  frequently  reported,  behavioral 
manipulations  can  result  in  changes  in  the  peak  latencies  of  visual 
primary  cortical  EP's  (3),  as  well  as  auditory  EP's  recorded  from  the 
round  window  (11).  The  anatomic  substrate  for  comparable  effects  in 
the  somatosensory  system  is  known  to  exist  (e.g. ,  26).  Thus,  EP 
changes  seen  after  Impact  may  be  assumed  to  be  related  to  injury  or 
precursors  of  injury  only  after  serious  consideration  of  these  other 
factors.  For  the  cervical  EP's  analyzed  here,  it  is  likely  that  there 
are  no  synapses  between  the  stimulating  and  cervical  recording  sites, 
and  that  the  observed  shifts  in  latency  result  from  changes  in 
conduction  velocities  of  spinal  fiber  groups.  Consequently, 
considerations  involving  cognitive  factors  should  not  apply.  However, 
since  small  increases  in  latency  often  occur  at  low  levels  of  impact, 
the  extent  to  which  the  observed  changes  are  within  the  limits  of  the 
"healthy"  operating  range  of  the  system  must  be  considered.  There  is 
no  simple  way  to  deal  with  this  problem,  and  therefore  the  approach 
taken  at  the  Naval  Biodynamics  Laboratory  is  to  evaluate  the 


*  Note  also  that  the  reported  temperature  relationships  would  require 
a  sudden  decrease  in  temperature  to  account  for  the  increase  in 
latency  resulting  from  Impact. 
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relationship  of  each  neurophysiological  measure  to  the  level  of  the 
impact.  Conclusions  concerning  CNS  dysfunction  depend  on  the  form  of 
the  relationship,  as  well  as  the  injury  threshold  determined  by 
clinical  and  neuropathological  examination. 

Three  broad  categories  of  changes  in  measures  of  CNS  function  may 
be  distinguished:  1  -  changes  observed  at  clearly  non-injurious  low 
levels  of  impact  (less  than  200  m/s2);  2  -  changes  that  occur  at 
intermediate  levels  of  impact  (200  -  800  m/s2);  3  -  changes  that  occur 
at  high  levels  of  probably  injurious  impact  (greater  than  800  m/s2). 
The  first  category  of  changes  may  be  excluded  from  further 
consideration.  This  includes  moderate  changes  in  the  amplitude  of  the 
cortical  slow  wave  and  shifts  in  the  latencies  of  the  cervical  EP 
peaks  less  than  25  jis/ms  (microseconds  shift  per  millisecond  peak 
latency).  The  latter  two  categories,  which  may  be  of  value  as 
neurophysiological  indices  of  injury  include:  shifts  in  the  latencies 
of  cervical  EP  peaks  exceeding  25  >is/ms;  large  changes  in  the 
amplitude  of  the  cervical  EP;  changes  in  the  shape  of  the  cervical  EP 
(2);  changes  in  the  cortical  EP  ripples;  and  substantial  and 
persistent  changes  in  the  cortical  primary  slow  wave  (19) .  In 
particular,  the  quantitative  analysis  of  shifts  in  latency  of  the 
cervical  EP  peaks  suggests  that  this  may  be  a  useful  measure.  The 
results  of  this  analysis  suggest  the  possibility  of  an  injury 
threshold  in  the  range  of  700  -  800  m/s2.  The  smaller  shifts  in 
latency  that  occur  near  600  m/s2  may  indicate  a  pre- injury  condition. 
This  is  an  encouraging  step  towards  identifying  criteria  for  assessing 
injury  to  the  CNS  caused  by  indirect  impact  of  the  head. 
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APPENDIX  A.  AEP's  From  All  Experiments  Analyzed 

Each  figure  In  this  appendix  shows  cortical  and  cervical  AEP's 
computed  with  N  -  50  from  a  single  experiment.  The  individual  figures 
are  referred  to  in  the  main  text  by  experiment  numbers,  each  of  which 
begin  with  the  characters  "LX"  (e.g. ,  the  first  experiment  number  is 
LX3008).  For  each  recording  site,  two  pre-impact  and  five  post-impact 
AEP's  are  arranged  in  a  column.  AEP's  in  the  same  row  are  taken  from 
the  same  time  epoch  relative  to  Impact.  The  time  indicated  in  each 
row  is  the  time,  relative  to  impact,  in  seconds,  of  the  stimulus 
resulting  in  the  first  EP  of  each  of  the  two  AEP's.  The  time  of  the 
stimulus  is  the  beginning  of  the  AEP.  Selected  constant  latency 
cursors  (vertical  axes)  are  shown  to  aid  in  the  observation  of  small 
shifts  in  latency.  The  number  at  the  top  of  each  cursor  is  the  median 
latency  (ms)  of  the  peak  from  AEP's  (N  ■  50)  computed  using  the  EP's 
from  the  two  minutes  Immediately  preceeding  impact. 
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APPENDIX  B.  Concepts  in  EP  Data  Analysis 
B-i.  Data  Space 

All  of  the  results  of  any  of  the  experiments  can  be  represented 
exactly  by  the  position  of  a  single  point  in  a  space  containing  a  very 
large  number  of  dimensions.  For  example,  if  the  spectral  properties 
of  the  EP  data  on  a  particular  channel  are  such  that  100  digital 
values  are  required  to  represent  each  EP,  and  there  are  1,500  EP's  on 
that  channel  in  the  experiment,  then  150,000  dimensions  would  be 
required  to  represent  that  portion  of  the  results.  This  is 
approximately  the  dimensionality  for  the  five  minutes  of  the 
post- impact  cervical  data  analysed  in  this  paper.  The  purpose  of  the 
data  analysis  is  to  reduce  the  dimensionality  to  a  small  number 
without  excessive  loss  of  the  information  that  is  crucial  to  the 
purpose  of  the  experiment. 

The  data  analysis  procedure  is  a  sequence  of  steps,  each  of  which 
reduces  the  dimensionality  by  one  or  more  orders  of  magnitude.  At 
each  step,  a  model  is  assumed,  which  if  correct,  will  insure  that  the 
information  loss  in  that  step  will  not  be  catastrophic.  The  success 
of  the  analysis  is  critically  dependent  on  the  applicability  of  the 
model.  For  example,  the  first  step  in  EP  work  is  usually  the 
computation  of  AEP's.  The  associated  model  includes  the  requirements 
that  the  important  information  is  contained  in  that  portion  of  the 
electrical  activity  that  is  synchronous  with  the  stimulus,  and  that 
the  time  course  of  the  important  changes  are  such  that  they  will  be 
observable  in  the  AEP’s.  Averaging  is  useful  only  to  the  extent  that 
these  are  true.  Averages  were  computed  with  N  »  10,  reducing  the  five 
minute  dimensionality  to  15,000. 

The  next  step  requires  the  reduction  of  the  dimensional ty  of  each 
AEP.  Though  the  means  for  accomplishing  this  remains  an  active  field 
of  research,  two  approaches  can  be  distinguished.  In  the  first 
approach,  the  extent  to  which  the  dimensionality  of  the  AEP  can  be 
reduced  without  loss  of  any  appreciable  amount  of  information  is 
determined.  In  the  second  approach,  information  other  than  EP  data  is 
available.  This  information  is  generally  in  the  form  of  different 
conditions  under  which  the  various  EP's  were  acquired.  In  this  case, 
the  approach  is  to  determine  the  extent  to  which  the  dimensionality  of 
the  AEP  can  be  reduced,  while  selectively  retaining  information  that 
is  related  to  the  different  conditions. 

In  the  present  work,  the  second  approach  is  taken.  Conditions 
are  assumed  to  be  changing  as  a  function  of  time  at  first  discretely 
at  the  moment  of  impact,  and  continuously  thereafter  as  a  consequence 
of  the  sequelae  to  Impact.  Various  algorithms  are  used  to  extract 
single  scalar  measures  from  each  AEP.  Scalar  measures  are  accepted  or 
rejected  on  the  basis  of  computational  stability  in  the  face  of  the 
variability  of  the  data,  as  well  as  the  ability  to  distinguish  the 
effects  of  Impact.  These  scalar  measures  Include  measures  of  the 
latency  and  amplitude  of  various  peaks,  and  quantitative 
representations  of  the  waveshape  of  various  portions  of  the  AEP.  Each 
scalar  measure  entails  its  own  set  of  assumptions.  In  general,  it  is 
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assumed  that  different  simply  measurable  aspects  of  the  AEP  will  be 
differently  related  to  Impact,  and  represent  different  functional 
effects  of  Impact.  In  the  present  report,  only  the  latency  shifts  of 
peaks  have  been  analysed  in  detail.  This  was  not  because  other  scalar 
measures  are  considered  to  lack  information,  but  because  of 
limitations  of  resources,  as  well  as  preliminary  indication  that  there 
was  Important  information  in  the  latency  shifts.  The  average  number 
of  peaks  measured  in  each  subject  was  about  six,  reducing  the  five 
minute  dimensionality  in  each  subject  to  900.  At  this  point,  each  of 
the  six  scalar  measures  (shifts  in  latency  at  the  six  peaks)  could  be 
plotted  as  a  function  of  time  allowing  observation  of  the  effects  of 
impact,  as  well  as  qualitative  comparisons  of  the  strength  and 
duration  of  the  effects  of  impact  among  subjects  and  conditions. 

These  observations  may  be  considered  to  be  subjective  reductions  of 
dimensionality.  For  more  objective  comparisons,  further  quantitative 
reduction  of  dimensionality  was  required. 


B-ii.  Exponential  Model  for  Reduction  of  Dimensionality 

To  reduce  the  dimensionality  further,  it  was  assumed  that 
following  Impact,  the  essential*  properties  of  the  time  course  of  each 
scalar  measure  could  be  represented  by  two  variables.  The  first  was 
the  amplitude  of  the  effect,  i.e. ,  a  measure  of  the  deviation  of  the 
the  scalar  measure  from  its  pre-impact  value  (not  to  be  confused  with 
the  amplitude  of  the  evoked  potential) .  The  second  was  the  duration 
of  the  effect.  At  this  point,  any  of  a  large  number  of  additional 
assumptions  could  be  made  to  allow  computation  of  the  variables.  It 
was  assumed  here  that  the  essential  effect  of  Impact  was 
instantaneous,  and  after  impact  it  decayed  exponentially  to  the 
pre- impact  level.  The  amplitudes  and  durations  were  usually 
determined  by  regression  of  a  single  exponential  decay  curve  on  each 
scalar  measure  from  the  first  five  post-impact  minutes.  In  most 
instances,  this  is  all  that  was  necessary.  In  some  cases,  however,  a 
substantial  portion  of  the  plot  of  the  scalar  measure  fell  below  the 
pre-impact  baseline  in  the  first  five  post-impact  minutes,  even  though 
the  initial  effect  was  an  Increase  in  the  scalar  measure.  In  these 
cases,  the  single  exponential  decay  model  was  grossly  inadequate,  and 
a  step  function,  with  the  step  occurlng  at  the  time  of  impact  was 
added.  Using  two-term  exponential  regression  the  effect  of  the  step 
was  minimized  and  a  single  amplitude  and  duration  could  still  be 
computed  in  most  cases. 

Assuming  six  peaks  per  subject,  the  five  minute  dimensionality 
for  each  subject  was  now  12.  The  dimensionality  was  finally  reduced 
to  two  by  computing  the  medians  of  the  amplitude  and  duration  across 


*  It  is  important  that  it  is  only  assumed  that  the  essential 
properties  of  the  scalar  measures  can  be  represented  by  the  two 
variables.  It  is  not  assumed  that  all,  or  even  most  of  the  details  of 
the  time  course  can  be  so  represented. 
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peaks  for  each  subject.  With  two  values  per  subject,  it  was  a  simple 
matter  to  construct  scatter  diagrams  for  observation  of  the  final 
experimental  result. 


B-iii.  Exponential  Model  for  Curve  Fitting 


Aside  from  determination  of  the  amplitude  and  duration,  it  was 
also  interesting  to  ask  whether  exponential  decay  would  be  a  good 
model  for  the  detailed  effects  of  Impact.  It  was  Immediately  apparent 
that  at  lower  Impact  levels,  the  measures  were  too  variable  to  provide 
an  answer.  Therefore,  only  data  from  high  level  impact  were  used  for 
this  purpose,  and  after  preliminary  work,  it  was  found  necessary  to 
consider  a  model  somewhat  more  complex  than  exponential  decay  alone. 

An  example  in  which  this  model  fits  the  data  reasonably  well  is  shown 
in  fig.  2.  The  complete  model  is  as  follows: 

1  -  Linear  baseline  shift.  It  was  assumed  that  in  the  absence  of 


Impact,  the  scalar  measure  may  be  changing  at  a  constant  rate 
with  time.  It  was  further  assumed  that  Impact  would  ..ot  alter 
the  process  responsible  for  this  change.  This  assumption  was 
based  on  the  observation  that  in  some  experiments,  some  of  the 
measures  appeared  to  change  in  an  approximately  linear  fashion 
during  the  experiment.  This  was  most  commonly  observed  with 
measures  of  latency,  which  tended  to  decrease  during  the 
experiment.  In  some  cases,  peak  latencies  tended  to  decrease 
during  an  entire  day  in  which  there  were  severed  experiments. 

For  example,  experiments  LX3008,  LX3009,  and  LX3010  were  all  run 
on  the  same  day  (table  1).  LX3008  and  LX3010  were  separated  by 

about  4  hours.  From  table  2  it  was  computed  that  the  average 
(across  the  four  peaks)  change  in  the  pre-impact  median 
latencies  was  -34  jis/ms  between  experiments  LX3008  and  LX3010. 

2  -  Impact-induced  shift  with  poly-exponential  decay.  It  was 

assumed  that  Impact  would  produce  an  instantaneous  shift  that 
decayed  back  to  baseline  according  to  a  poly-exponential 
function  decay  (the  sum  of  several  simple  exponential  decay 
terms).  The  simplest  assumption  for  a  process  that  appears  to 
be  decaying  with  time  toward  an  asymptote  is  that  the  decay  is  a 
single  exponential  one.*  Some  of  the  plots  of  the  scalar 
measures  as  a  function  of  time  suggested  this  possibility. 


*  This  amounts  to  assuming  that  the  rate  of  decay  at  any  given  time 
(relative  to  impact)  is  proportional  to  the  remaining  strength 
(amplitude)  of  the  effect  at  that  time.  The  constant  of 
proportionality  is  the  reciprocal  of  the  time  coefficient  of  the 
decay.  This  simple  description  is  applicable  to  a  large  number  of 
physical,  chemical,  and  physiological  processes.  For  example.  Brown 
and  Brown  (4)  report  that  changes  in  cardiovascular  measurements 
produced  by  a  direct  occlputal  blow  in  Rhesus  returned  to  normal  with 
an  exponential  time  course. 
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However,  In  some  cases  it  appeared  that  partial  recovery  took 
place  quite  rapidly,  while  part  of  the  effect  of  Impact  took  a 
much  longer  time  to  recover,  suggesting  the  poly-exponential 
model. 

3  -  Impact-induced  "permanent”  shift.  It  was  assumed  that  the 

Impact  would  produce  an  instantaneous  change  in  the  measure 
persisting  unchanged  for  the  duration  of  the  experiment.  In 
practice,  it  was  found  that  in  the  least  noisy  data,  once  the 
magnitude  of  a  time  coefficient  exceeded  several  times  the 
duration  of  the  interval  analyzed,  it  became  impossible  to 
distinguish  a  slow  decay  from  a  permanent  shift. 

4  -  Residual  noise.  This  is  the  variation  not  accounted  for  by  the 

above  three  parts. 

The  model  can  be  represented  by  the  following  equation: 

N 

y  -  St  +  B  +  h(t)D  h(t)aiexp(t/Ti)  +  E(t) 

i-1 

where: 

y  is  the  value  of  the  measure, 

t  is  time  relative  to  impact, 

S  is  the  slope  of  the  linear  baseline  shift, 

B  is  the  amplitude  value  of  the  baseline  Instantaneously 
prior  to  Impact, 

D  is  the  "permanent"  shift  induced  by  Impact, 

N  is  the  number  of  exponential  terms, 

a^  are  the  exponential  amplitudes, 

T^  are  the  exponential  time  coefficients, 

h(t)  is  the  unit  step  function  with  value  zero  for  t  less  than 
zero,  and  unity  for  t  greater  than  or  equal  to  zero, 

E(t)  is  the  noise. 


B-iv.  Regression  Procedures  for  Curve  Fitting 

The  linear  baseline  shift  parameters,  S  and  B,  were  determined  by 
applying  simple  linear  regression  to  the  pre-impact  data.  The  linear 
equation  was  extrapolated  into  the  post- Impact  region,  and  values 
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computed  from  the  linear  equation  were  subtracted  from  the  data.  The 
residual  was  then  subjected  to  poly-exponential  regression,  providing 
estimates  of  the  a^  and  T^. 

The  poly-exponential  regression  was  also  applied  to  the  residual 
to  estimate  the  "permanent"  shift  due  to  Impact,  D.  This  was  done  by 
using  an  additional  exponential  term  in  the  regression  procedure.  If 
there  were  a  "permanent"  change  in  the  measure,  the  extra  exponential 
term  would  have  a  time  coefficient  (T)  of  either  sign  that  was  very 
long  in  relation  to  the  duration  of  the  analysis  interval.  The  sign 
of  the  amplitude  of  this  term  depended  on  the  direction  of  the 
permanent  shift.  For  all  values  of  t  during  the  analysis  interval, 
t/T  is  close  to  zero  and  exp(t/T)  is  close  to  unity  irrespective  of 
the  sign  of  T.  The  corresponding  amplitude  is  thus  an  estimate  of  D, 
the  "permanent"  shift  induced  by  impact.  In  practice,  the  regression 
algorithm  did  not  always  produce  a  term  with  a  long  time  coefficient, 
especially  if  the  "permanent"  displacement  was  small  or  non-exlstant, 
or  if  the  "permanent"  shift  was  in  the  same  direction  as  the 
exponential  decay.  Rather,  it  used  the  additional  term  to  fit  some 
detail  of  curvature,  triiich  usually  resulted  in  substantial  distortion 
of  the  other  coefficients.  However,  the  algorithm  could  often  be 
induced  to  provide  the  term  with  a  long  time  coefficient  by  simply 
adding  a  constant  value  to  all  of  the  data.  The  permanent  shift,  D, 
was  then  the  difference  between  the  added  constant  and  the  amplitude 
of  the  exponential  term  with  the  long  time  coefficient.  If,  in  fact, 
there  were  no  "permanent"  shift  at  Impact,  this  difference  would  be 
very  small. 


B-v.  The  Distribution  of  the  Latency  Shifts  —  Outliers 

Generally,  tfien  a  particular  type  of  result  is  to  be  derived  from 
a  given  data  set,  some  identifiable  subsets  of  the  available 
information  are  more  reliable  than  others.  Therefore,  various 
statistical  procedures  have  been  developed  which  selectively  use 
information  that  has  a  high  probability  of  being  reliable.  A  common 
example  is  the  use  of  the  median  to  estimate  central  tendency.  In 
this  case,  only  part  of  the  available  Information  is  used;  the  actual 
magnitudes  of  most  of  the  data  values  are  discarded.  The  elimination 
of  outliers  provides  another  example.  In  the  current  analysis,  both 
of  these  procedures  are  combined,  resulting  in  a  substantial 
improvement  of  the  reliability  of  the  estimates  of  the  amplitude  and 
duration  of  the  observed  effects. 

Poly-exponential  regression  is  very  sensitive  to  outliers.  This 
sensitivity  arises  from  two  sources.  First,  the  squaring  in  the  mean 
square  deviation  criterion  amplifies  the  affect  of  outliers.  Second, 
there  are  strong  correlations  among  the  regression  coefficients  which 
allow  relatively  small  changes  in  the  data  to  affect  the  coefficients 
sufficiently  so  that  the  essential  character  of  the  result  is  changed. 
This  often  happens  without  an  appreciable  change  in  the  mean  square 
deviation.  The  problem  is  especially  acute  if  more  than  one 
exponential  term  is  used.  Using  poly-exponential  regression,  however, 
it  is  possible  to  devise  a  procedure  that  eliminates  outliers. 
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Before  selecting  outliers,  some  understanding  of  the  distribution 
of  the  data  Is  necessary.  Assume,  for  the  moment,  that  there  Is  no 
Impact,  and  that  the  time  series  of  evoked  potential  vectors  Is 
stationary.  In  the  case  of  latency  measurements  (either  by  simple 
peak  detection  or  normalized  cross-correlation)  the  distribution  of 
the  measurements  may  be  considered  as  the  sum  of  three  distributions 
as  follows: 

1  -  Correct  detection  of  the  target  peak.  The  distributions  of  the 

target  latency  measurements  are  assumed  to  resemble  normal 
distributions.  In  the  absence  of  experimental  manipulation,  in 
all  cases  \rtiere  the  cervical  EP  Is  clearly  present,  it  was 
found  that  the  latencies  for  the  various  peaks  are  very  stable 
over  short  periods  of  time.  The  standard  deviations  of  the 
latency  distributions  are  of  the  order  of  50  microseconds  or 
less. 

2  -  Absence  of  the  Evoked  Potential .  Due  either  to  artifact,  or  to 

true  physiological  absence  of  the  EP,  the  expected  distribution 
of  the  latency  measurements  is  rectangular  over  the  latency 
range  of  the  peak  search.  Since  the  peak  is  generally  not 
found  in  the  middle  of  the  range,  this  rectangular  distribution 
is  asymmetric  with  respect  to  the  mean  of  the  distribution  of 
target  detections.  A  typical  width  of  the  search  range,  and 
thus  the  width  of  the  rectangular  distribution,  is  500  to 

1000  ps . 

3  -  Detection  of  an  incorrect  peak.  This  distribution  will  depend 

on  the  latencies  of  available  incorrect  peaks.  Often,  there  is 
one  extra  peak  in  the  search  interval.  The  resulting 
distribution  is  severely  asymmetric  with  respect  to  the 
theoretical  mean  of  the  correct  target  distribution.  Rarely, 
two  extra  peaks  are  present.  The  distribution  still  results  in 
asymmetry,  because  the  peaks  do  not  occur  at  regular  intervals, 
and  are  not  detected  with  equal  probability.  The  smallest 
Intervals  between  peaks  in  the  cervical  data  are  of  the  order 
of  300  ps. 

Thus,  the  observed  distribution  has  excessively  large  and  asymmetric 
tails,  as  compared  with  the  target  distribution.  It  is  clear  that 
because  of  the  large,  asymmetric  contributions  of  the  two  "error" 
distributions,  the  mean  of  the  observed  distribution  will  not  be  a 
good  estimate  of  the  target  mean.  The  median  will  provide  a 
substantially  better  estimate  of  the  target  mean  if  more  than  50%  of 
the  peak  detections  are  target  peak  detections. 

Because  of  the  non-statlonarity  due  to  impact,  the  above 
considerations  cannot  be  directly  applied,  but  the  following  procedure 
may  be  used  to  eliminate  outliers  from  the  post-impact  data  if  the 
number  of  outliers  is  not  too  great: 

1  -  Compute  a  least-squares  poly-exponential  regression  using  the 
entire  data  set. 
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2  -  Subtract  the  poly-*exponential  function  from  the  data  leaving  a 

"residual". 

3  -  Discard  outliers  relative  to  the  median  of  the  residual  using 

some  predetermined  criterion. 

4  -  Recompute  the  poly-exponential  regression  using  the  remaining 

data. 

In  this  procedure,  it  is  assumed  that  the  post- Impact  data  may  be 
considered  to  be  the  sum  of  a  continuation  of  the  stationary 
pre-impact  time  series  and  the  effect  of  Impact  on  the  time  series. 

The  short  term  stability  of  the  data  are  assumed  to  be  unchanged  by 
the  Impact  except  for  this  additive  effect.  It  is  further  assumed 
that  the  initial  poly-exponential  regression  is  close  enough  to  the 
actual  effect  of  impact  so  that  the  data  in  the  tails  of  the  residual 
distribution  can  be  considered  outliers.  Since  the  "error" 
distributions  are  much  wider  than  the  target  distribution,  a  very 
rough  initial  fit  will  accomplish  this,  as  was  judged  by  examination 
of  graphs  of  the  initial  regression  and  the  associated  residual . 

Procedures  for  objective  elimination  of  outliers  are  in  the 
planning  stage.  In  the  current  analysis,  outliers  were  removed 
subjectively  after  exhaustive  examination  of  the  graphic  results  of 
initial  poly-exponential  regressions.  Sometimes  the  procedure  was 
iterated  several  times.  The  final  results  (tables  2  &  3,  figs.  4  &  5) 
were  obtained  using  the  following  guidelines: 

1  -  Any  coefficient  resulting  from  a  few  data  points  was 

eliminated,  generally  by  discarding  the  data  points.  Thus  no 
time  coefficients  with  magnitudes  less  than  10  seconds  are 
reported. 

2  -  Where  the  data  were  clearly  multimodal,  all  data  outside  of  the 

main  mode  were  discarded. 

3  -  Where  the  data  were  highly  variable,  but  not  clearly 

multimodal,  no  coefficients  were  computed. 

4  -  Where  shifts  in  latency  were  small  and  no  apparent  pattern  was 

seen,  various  procedures  were  used  to  produce  a  regression  that 
would  represent  most  of  the  plot.  The  time  coefficients  in 
such  cases  are  meaningless. 

5  -  Positive  time  coefficients  indicating  exponential  growth  were 

considered  meaningless  in  the  current  context.  Unless  they 
were  so  large  that  the  contribution  of  the  term  was  essentially 
constant  over  the  time  Interval  under  study,  they  were 
eliminated  by  various  means. 

6  -  The  first  few  points  after  impact  were  eliminated.  The 

exponential  amplitudes  therefore  represent  values  that  are 
extrapolated  back  from  post- Impact  data. 
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7  -  Last,  but  possibly  most  significant,  data  were  discarded  if 
they  indicated  shifts  in  latency  larger  than  any  seen  in 
clearly  resolved  AEP’s  under  comperable  conditions.  This  was 
especially  Important  in  high  level  impact  experiments.  It  is 
considered  that  there  is  some  uncertainty  regarding  the 
magnitude  of  the  largest  shifts  in  latency  occurlng  under  these 
conditions. 


B-vi.  Measurement  of  Duration  —  Comparison  With  Other  Methods 

In  past  work  the  duration  has  been  emphasized,  using  "threshold 
crossing"  procedures  (19,  32).  These  measure  the  time  between  Impact 
and  the  first  value  of  the  scalar  measure  to  return  to  some  threshold 
level,  such  as  the  pre-impact  median,  or  one  standard  deviation  from 
the  pre-impact  grand  mean.  With  moderately  noisy  data  such  as  ours, 
this  simple  procedure  may  be  misleading  unless  some  sort  of  smoothing 
is  applied  to  the  data  before  measuring  the  time  to  threshold 
crossing.  However,  if  the  noise  is  of  such  a  nature  that  outliers  are 
frequent,  the  smoothing  procedure  can  be  further  improved  by 
incorporating  a  method  that  selectively  discards  outliers.  The 
exponential  regression  method,  as  used  here,  is  a  smoothing  procedure 
of  this  type. 

Any  procedure  for  discarding  outliers  during  the  non-stationary 
period  following  Impact  must  involve  assumptions  about  the  effect  of 
Impact.  The  exponential  model  Implies  that  except  at  the  instant  of 
impact,  there  are  no  sudden  changes.  Other  models  may  be  used  to 
exclude  outliers,  and  it  is  likely  that  many  would  do  at  least  as  well 
as  the  exponential  decay  model.  The  final  results  should  be 
relatively  insensitive  to  the  choice  of  the  model  for  the  same  reason 
that  the  exponential  model  works  even  if  the  exponential  coefficients 
are  known  only  approximately:  the  error  distributions  are  much  wider 
than  the  target  distributions. 


